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Abstract. - 

A theoretical formula for the probability density function (PDF) of velocity derivatives in a 
fully developed turbulent flow is derived with the multifractal aspect based on the generalized 
measures of entropy, i.e., the extensive Renyi entropy or the non-extensive Tsallis entropy, and 
is used, successfully, to analyze the PDF's observed in the direct numerical simulation (DNS) 
conducted by Gotoh et al.. The minimum length scale rd/n in the longitudinal (transverse) 
inertial range of the DNS is estimated to be r\/n = 1.716 (rj/n — 2.180) in the unit of the 
Kolmogorov scale n. 



In the previous paper 0, we analyzed, precisely, the probability density functions (PDF's) 
of velocity fluctuations observed in the direct numerical simulation (DNS) of turbulence con- 
ducted by Gotoh et al. j^] at the Taylor microscale Reynolds number R\ = 381. At this 
Reynolds number, the PDF's had been measured with high accuracy up to the order of 
10 -9 ~ 10~ 10 . The correctness of the analytical formula for the scaling exponents of velocity 
structure function |^-^| enabled us to extract the value of the intermittency exponent \x by 
fitting it with the ten observed data in the DNS by the method of least square 0. With 
the intermittency exponent, the parameters in the analytical formula of the PDF of veloc- 
ity fluctuations are determined, self-consistently. The formula was used to fit each of 
the observed PDF's of velocity fluctuations for ten different separations r/nhy means of the 
method of least square [JjJ, and to extract, successfully, the dependence of the number n of 
steps in the energy cascade on the separation r/iy. 

n = -1.050 x log 2 r/rj + 16.74 (for % < r), (1) 
n = -2.540 x log 2 r/n + 25.08 (for r < l\) (2) 

with the crossover length i^/rj = 48.26 for longitudinal fluctuations, and 

n = -0.9896 x log 2 r/rj + 13.95 (for ij < r), (3) 
n = -2.820 x log 2 r/n + 23.87 (for r < ^F) (4) 
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with the crossover length tj/rj = 42.57 for transverse fluctuations These straight lines on 
semi-logarithmic sheet, explicitly, told us that there exist two scaling regions i.e., the upper 
scaling region with larger separations which may correspond to the scaling range observed by 
Gotoh et al. and the lower scaling region with smaller separations which is another scaling 
region extracted first by the systematic analyses in PJ. These scaling regions are divided 
by the crossover lengths &/r) and jr\ approximately of the order of the Taylor microscale 
A/17 = 38.33 reported in || at R\ = 381. 

We can view that the turbulent flow, satisfying the Navier-Stokes equation 

du/dt+{u-V)u = -V {p/p) + vV 2 u (5) 

of an incompressible fluid, consists of a cascade of eddies with different diameter of the order 
of £ n = 5 n £o where S n = 2~ n {n = 0, 1, 2, • • •). The quantities p and p represent, respectively, 
the mass density and the pressure. At each step of the cascade, say at the nth step, mother- 
eddies of size £ n -i produce smaller daugter-eddies having a half of mother's diameter with the 
energy-transfer rate e„ that represents the rate of transfer of energy per unit mass from eddies 
of size l n -\ to those of size £ n (the energy cascade model). Following formally the energy 
cascade model, we are measuring space scale by l n . However our analysis in the following is 
not restricted within the energy cascade model, i.e., the number of step n can be real number. 

In this paper, we will derive the formula for the PDF of velocity derivatives in fully_ de- 
veloped turbulence by the statistics based on the generalized entropy, i.e., Renyi's 




Tsallis' 11, 12j, and will analyze the velocity derivative PDF's obtained in the DNS 
R\ — 381 having far better accuracy than that in any previous experiments, real or nu- 
merical. The main interest, here, is the longitudinal velocity derivative du{r) / dr\ and the 
transverse velocity derivative du{f)/dr^ (or du{r)/dr^,), where r = (ri,r2,rs), and u is the 
?*i -component of the fluid velocity field u of the turbulent flow produced by a grid with size 
£0 putting in a laminar flow parallel to the r*i direction. Introducing the velocity difference 
Su n of the component u at two points separated by the distance £ n , the velocity derivatives 
may be estimated by \\mi n ^ 8u n / £ n (= linin^oo 5u n /£ n ). The Reynolds number Re of the 
system is given by Re = 5uq£q/v — (£o/rj) 4 ^ 3 . For high Reynolds number Re 3> 1, or for the 
situation where effects of the kinematic viscosity v can be neglected compared with those of 
the turbulent viscosity, the Navier-Stokes equation (|5|) is invariant under the scale transfor- 
mation Jl|,[l!l: f -> Ar, u -> A a/3 u, t -> A 1 " Q / 3 t and {p/p) -> A 2q / 3 (p/p). The exponent a 
is an arbitrary real quantity which specifies the degree of singularity in the velocity deriva- 
tive JT^ I for a < 3, i.e., lim^^o 5u n /£ n ~ lim^ n _,o GJ Z 1 which can be seen with the relation 
5u n /8u = (V4) Q/3 - 

The energy E n = (Su n ) 2 /2 per unit mass contained in an eddy of size £ n is estimated 
as E n ~ (£ n uj n ) 2 where w„ ~ Su n /£ n represents the angular momentum of the eddy. We 
see that lim^^o^n has the same singularity as the velocity derivative does. Within the 
region satisfying the scale invariance, we have E n — EoS n a ^ 3 . The energy spectrum E(k), 



defined through E n = J^ n+1 dkE{k) with k n = £ n 1 , has the wavenumber dependence E(k) cx 
£,-i-2q/3 gy wa ^ energy transfer rate, estimated by e„ ~ E n uj n , satisfies e„ = eo5"~ x 
in the region of the scale invariance. Kolmogorov's assumption in K41 |J that there is no 
fluctuation in e n leads us to a = 1. With this value of a, E(k) represents the Kolmogorov 
energy spectrum, i.e., E(k) oc fc~ 3 / 5 . In order to explain the intermittency in turbulence, we 
shall introduce a fluctuation in a. 

( 1 ) Here, rj is the Kolmogorov scale |J defined by rj = (i^/e) 1 / 4 with u being the kinematic viscosity, and 
t (= en) the energy input rate to the largest eddies with size £q, and has the value r\ = 2.58 X in the 

DNS B| at R\ = 381. 
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The present analysis rests on the assumption that the distribution of the exponent a is 
multifractal, and that the probability PW {a)da to find, at a point in physical space, an eddy 
of size £ n having a value of the degree of singularity in the range a ~ a + da is given by 
p(»)(a) = [pW(a)] n with 

pW(a)oc [l-( a -a ) 2 /(4a) 2 ] 1/(1 - 9) , (Z\a) 2 = 2A/[(1 - g) In 2]. (6) 
Here, it is assumed that each step in the cascade is statistically independent. The distri- 



bution function (g) is derived by taking an extremum of the generalized entropy Q with 
the two constraints, i.e., the normalization of distribution function: J daP^(a) = const, 
and the q- variance being kept constant as a known quantity: er 2 , = (J daP^ 1 '(a) q (a — 
a ) 2 )/ JdaPW{a)i. 

The dependence of the parameters ao, X and q on the intermittency exponent fi is de- 
termined, self-consistently, with the help of the three independent equations, i.e., the energy 
conservation: (e„) = e, the definition of the intermittency exponent /i: (e 2 ) = e 2 5~^, and 
the scaling relation K 3 "!: 1/(1 — q) = l/a_ — l/a + with a± satisfying f(a±) — where the 
multifractal spectrum || [J 

/(a) = 1 + (1 - qT 1 log 2 [l - (a - a f / (Aa) 2 ] (7) 

is derived by the relation P( n \a) oc Sn |], 13 that reveals how densely each singularity, 
labeled by a, fills physical space. The average (■ ■ •) is taken with P^(a). Note that the 
relation between a and e„ is given by e„/e = (£ n /£o) a ~ ■ For the region where the value of 
fi is usually observed, i.e., 0.13 < [i < 0.40, the three self-consistent equations are solved to 
give the approximate equations Q]: a = 0.9989 + 0.5814^, X = -2.848 x 10 -3 + 1.198/x and 
q = -1.507 + 20.58^ - 97.11/i 2 + 260.4^ 3 - 365.4^ 4 + 208.3^ 5 . 

Let us suppose that id is the typical length giving the minimum scale within the energy 
cascade model. With this shortest length ^d, the velocity derivatives may be given by \s\ = 
Su^/id with Sud = Su n=rid where rid is introduced through £d/£a = $d = 2~" d . Since there are 
two mechanisms in turbulent flow to rule its dissipativc evolution, i.e., the one controlled by 
the kinematic viscosity that takes care thermal fluctuations, and the other by the turbulent 
viscosity that is responsible for intermittent fluctuations related to the singularities in velocity 
derivative, it may be reasonable to assume that the probability lJ(x)dx to find the scaled 
velocity derivative x = Sdtos in the range x ~ x + dx can be divided into two parts: 

n(x)dx = n s (\x\)dx + An(x)dx. (8) 

Here, the singular part PDF _/7g(|a;|) represents the contribution from multifractal distribution 
of the singularities, and the correction part AII(x) from the viscus term neglected in the scale 
transformation. The former is derived through IIs(\x\)dx = pw (a)da with the transforma- 
tion of the variables: \x\ = S^ 3 ■ The mth moments of the velocity derivatives, defined by 



x 



= J_ oo dx\x\ m II (x) , are given by 



= 2 7m + (1 - 2 7o ) a m S C d m (9) 



( 2 ) The Renyi entropy S^pW (a)] = (1 - q)' 1 In J daPW(a)« |u| has the extensive character as the usual 
thcrmodynamical entropy does, whereas the Tsallis entropy [P^ (a)] = (1 — (J" da (a) 9 — 

12l 19] is non-extensive. In spite of different characteristics of these entropies the distribution functions giving 
tne extremum of each entropy have the common structure (n). 

( 3 ) The scaling relation is a generalization of the one derived first in Jl7| , |l8) to the case where the multifractal 
spectrum has negative values. 
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with a 3q = {2/[C\ ,2 {l + C g - /2 )]}i/2 7 Cq = l + 2f{l - q)X In 2 and 2 7m = dx \x\ m II N {x). 
We used the normalization: ((1)) = 1. The quantity 

Cm = ® m/3 — 2Xm 2 / 

is the so-called scaling exponent of the velocity structure function, whose expression was 
derived first by the present authors ||-[| . Note that the formula is independent of the length 
£d, and, therefore, independent of £ n . 



longitudinal transverse 





0.240 


0.327 


q 


0.391 


0.543 


ato 


1.138 


1.189 


X 


0.285 


0.388 



Table I - Values of the intermittency exponent /i and the parameters q, eta and X for longitudinal and 
transverse velocity fluctuations |lj, determined by the formula being consistent with the observed 
data in DNS conducted by Gotoh et al. 



9 1 + C 



1/2 
m/3 



^(i + c^) /(i -e) 



(10) 



With the help of the analytical formula ([H]), we determined in Q| the values of the in- 
termittency exponent /x and the parameters g, «o & n d X by fitting the ten DNS data of the 
scaling exponents Cm ( m — T 2, • • • , 10) at Ra = 381 with the method of least squares. 
The determined values are listed in table | both for the longitudinal and transverse velocity 
fluctuations. Note that the relation /i = 2 — ^ is satisfied within the experimental error bars. 
We have a+ - a a = a - a- = 0.6818 (0.8167), Aa = 1.160 (1.566) for the longitudinal 
(transverse) fluctuations. 

Since we are interested in the large deviation stemmed from the singular part 77s(|a;|) that 
may contribute to the symmetric part of the PDF, we will symmetrize right and left of the 
experimental PDF, and will compare it with our theoretical PDF in the following. 

Let us introduce the PDF /7(ICI) °f the velocity derivatives by /7(|£|)cf£ = II(\x\)dx with 
the new variable 

z = s /({s*y/ 2 = x /((x 2 )}^, iei = ^ d Q/3 - C2/2 , (ii) 

scaled by the variance of velocity derivatives. It may be appropriate to devide the PDF into 
two parts: 

n(0 = iz*<(0 for e < |c| < ^™ /3 - C2/2 (12) 

# (0 - n<*(0 for ici < r (13) 

with the point £* defined by = ^*/ 3 " C2 / 2 where a* is the solution of C 2 /2-a/3+l -/(a) = 
that provides us with the least nd-dependence of Here, £ = [272^+ (1— 2j )a2]~ 1 / 2 ■ 

As the value of £* turns out to be of order 1 (see below), we are deviding the region by the 
order of the variance of velocity derivative. #<*(£) and /?*<(£) are connected at £* under the 
condition that they should have the same value and the same slope there. In our analysis, since 
we are assuming that the large deviations of the velocity derivative come from the maltifractal 
distribution of these singularities in real space |Xq] , it is consistent to put //*<(£) = Z7s(|£|). 
This leads 

- r 9 *i n d/ (1— q) 

n*<(0 = n s (e/|e|) {l - [31n|^ | / (Aa | ln<5 d |)] 2 } (14) 
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with |£o| = ^d°^ 3 ^ 2 ^ 2 - For smaller velocity derivatives |£| < £*, the contribution to the 
PDF may come, mainly, from the term ignored in the scale transformation producing, i.e., 
the term responsible for thermal fluctuations related to the kinematic viscosity, and also from 
measuremant errors, we take for the PDF /7<*(£) a Gaussian function. The connection at £* 
with (jiil) gives us 

£<*(0 = ffge-W^HK/n'-il/s (15) 

with 5s = 3(1 - 2j )/(2^2ttX\ ln<5 d |). It is remarkable that the PDF Z7(|f|) of the velocity 
derivative, given by ( fL2[ ) and p3; ) with (14) and (|l5), turns out to have the same structure as 
the PDF ij( n )(|£„|) of the velocity fluctuations jl|,|6-|| with the separation ^„ = 




Fig. 1 - PDF's of the velocity derivatives for (a) longitudinal and for (b) transverse. Closed circles 
are the symmetrized points of the PDF observed in DNS conducted by Gotoh et al. at R\ — 381 ||. 
Solid line represents the PDF given by the analytic formula (^) and ( |l3| ) with the parameters given 
in Table and with (a) rid = = 23.1 and (b) = nj = 20.7. 

The comparison between the PDF's of the velocity derivatives at R\ = 381 measured in 
the DNS Q and those obtained by the present analysis is given in fig. [l] (a) for longitudinal 
derivative and in fig. [I] (b) for transverse derivative. In order to extract the symmetrical 
part of the PDF, we took mean average between the DNS data on the left hand side and 
that on the right hand side. The symmetrized data are described by closed circles. The 
solid lines are the curves of /?(£) given by ( |l4| ) and ( fi"5| ) with the values of parameters in 
table H Note that £* = 0.982 for the longitudinal velocity derivative, whereas £* = 0.900 
for the transverse derivative. The number — 23.1 (nj — 20.7) of steps in the cascade for 
longitudinal (transverse) velocity fluctuations is derived by the method of least squares with 
respect to the logarithm of PDF's for the best fit of our theoretical formulae, consisting of ( |l4|) 
and ([15]), to the observed values of the PDF by discarding those points which have observed 
values less than 10~ 10 since they scatter largely in the logarithmic scale. We see an excellent 
agreement between the measured PDF for the velocity derivatives and the analytical formula 
of PDF derived by the present self-consistent theory. 

Substituting the obtained values n\ = 23.1 (nj — 20.7) of the number of steps in the 
energy cascade for longitudinal (transverse) velocity fluctuations into (0) (into we obtain 
the shortest length rd of separation in the inertial range with the value r\jr\ = 1.716 {rj /r/ = 
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2.180). Then, we can conclude that the range of the lower scaling region in the inertial 
range of the longitudinal (transverse) velocity fluctuations is given by rWr/ < r/rj < l\jr\ 
( r d fv — r / f l — ^c/v)- Adding to this the upper scaling region, the total inertial range 
for longitudinal (transverse) fluctuations within the DNS Q at R\ — 381 turns out to be 
1.716 < r/rj < 1220 (2.180 < r/r, < 1220), where the value 1220 is the largest separation 
taken by Gotoh et al. Q for the measurement of the PDF's of velocity fluctuations. Note 
that the largest scale £0 for the longitudinal (transverse) fluctuations can be estimated by 
putting n = into (0) (into (§)), i.e., = 6.299 x 10 4 (J%/r) = 1.752 x 10 4 ) leading to 

Rc L = 2.506 x 10 6 (Rer = 4.550 x 10 5 ). We see that the shortest length scale £ d /r, = 2~™ d V? 
for longitudinal (transverse) velocity fluctuations within the energy cascade model becomes 
l\jr] = 7.006 x 10~ 3 (£j/rj = 1.029 x 10" 2 ) which is different from the length r^/rj (rj /if) 
giving the estimate of the lowest end of inertial range in the DNS Q at R\ = 381. We 
presume, here, that r/rj in the formulae (|l|), (|2|) and (||), (^) provides us with a real distance 
in the support of the velocity fields u(r). 

Summarizing, we derived in this paper the formula for PDF of velocity derivatives, (O) 
and (13) with ( |l4| ) and (|l5|), and showed that it explains the observed PDF's in the DNS J^, 
precisely up to the order of 10" 10 , with the parameters in Table 0. The latter analysis provides 
us with the number rid of steps in the energy cascade. The shortest length scale (rj), which 
serves the lowest end of the inertial range, is derived by making use of the obtained value of 
n d ( n d ) by assuming that the formula (g) (the formula (jij)) between n and r/rj, derived in the 
analyses of the PDF of longitudinal (transverse) velocity fluctuations jj], is applicable even 
for this shorter scale of length. 

Let us close this paper by mentioning something about another trial for deriving the PDF 
of velocity derivatives. We introduced, as a preliminary test, the same philosophy for the 
definition of velocity derivative proposed by Benzi et al. |ll| into the present analysis based 
on the generalized statistics, and saw that the PDF of velocity derivatives thus obtained (see 
(43) in Q) cannot explain the PDF's provided in the DNS ||] so precise as the present PDF 
given by (12) and (fill). A further investigation on the relation between the approach in |ll| 
and the present one may be one of the interesting future problems. 

The authors are grateful to Prof. T. Gotoh for enlightening discussion and his kindness to 
show his data prior to publication. The authors would like to thank Prof. C. Tsallis for his 
fruitful comments with encouragement. 
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